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The paper discusses the single-mode Jaynes-Cummings model with time dependent parameters. 
Solvable models for two-level systems are utilized to consider the changes in the photon distribution 
affected by the passage of atoms through the cavity, ft is suggested that such systems may be 
used as filters to modify the photon distribution. The effect can be enhanced by repeatedly sending 
new atoms through the cavity. We show that such filters can cut out either small or large photon 
numbers, ft is also shown that the method can be used to narrow down photon distributions and 
in this way achieve highly non-classical sub-Poissonian states. Some limitations and applications of 
the method are presented. 

PACS numbers: 



I. INTRODUCTION 

Cavity QED is based on the assumption that the mode 
selection allows us to treat the radiation modes as a dis- 
crete set and consider states with only a few modes ex- 
cited. Thus the system is described by a finite number 
of degrees of freedom; the rotating wave approximation 
(RWA) sets further limits to the states which have to be 
included. In the case involving a single mode only, we find 
the solvable Jaynes-Cummings model, which has served 
as a prototype for many physical situations in Quantum 
Optics. 

The Jaynes-Cummings model was introduced in 1963 
1] and it has since played a central role in the develop- 
ment of Quantum Optics; see the review 0. It has also 
served as the theoretical foundation for analysing the re- 
sults of experiments in cavity QED Q . 

In all these applications, however, the model has been 
treated with constant parameters. In our earlier work 
we discussed the model when the parameters are allowed 
to change gradually during the interaction. We now 
propose to continue our work on such time dependent 
phenomena in the framework of the Jaynes-Cummings 
model. 

In this paper, we assume that we can consider such 
a case where two-level atoms consecutively interact with 
a single cavity mode sharp in frequency. The various 
decay constants are taken to be zero. Thus all time scales 
involved must be shorter than the damping times of the 
physical degrees of freedom. 

In Sec. m we introduce the basic Jaynes-Cummings 
model with time dependent coefficients. We also discuss 
the possible solutions and the validity conditions for the 
approach. In Sec. IHII we illustrate the filtering action of 
atom dynamics by the simple results following from the 
Landau-Zener model. The actual discussion of the pho- 
ton state filtering is given in Sec. IIVI where we apply the 
solution of the more realistic Demkov-Kunike model. We 
find that this may be utilised both to achieve sharpening 
of the photon distribution and, if we prefer, as a low or 



high pass photon number filter. With numerical illus- 
trations, we show that both narrowing of the state and 
low pass photon filtering may be achieved. In the former 
case, we show that highly non-classical, sub-Poissonian 
photon distributions can be obtained. We also present 
some approximate analytical expressions for the width of 
the ensuing distributions. These are shown to provide 
god approximations to the numerically obtained exact 
ones. Finally we conclude the paper in Sec. with a 
summary and some possible uses of these methods. 



II. ADIABATIC JAYNES-CUMMINGS MODEL 

In this paper, we consider the Jaynes-Cummings model 
with time dependent parameters. There are two essential 
ones; within the RWA they are the cavity detuning and 
the coupling constant. Both are to a large extent un- 
der experimental control in modern experiments. Here 
we consider the situation when these are made to change 
slowly during the course of the interaction between the 
cavity and atoms. One assumption here is, of course, that 
the timing of the atoms in the cavity can be controlled 
by the experiment. This is non-trivial but not unreal- 
istic. The change must be slow, because we want the 
quantum state to follow the change adiabatically. This 
is necessary in order to preserve the identity of the single 
radiation mode involved; too fast a change will mix in 
higher modes and destroy the simplicity of the system. 
Thus we are working with an adiabatic extension of the 
ordinary cavity QED system. 

The Jaynes-Cummings Hamiltonian separates into 
two-component families of coupled states, and we assume 
that this separation holds even when the parameters are 
taken to vary slowly enough. Each family is characterised 
by an initial photon number in the cavity, and the total 
solution emerges from a sum over all initial photon states. 
The problem now becomes one of coupled two-level sys- 
tems with time dependent coefficients. Choosing these 
in suitable ways, we can effect a variety of transforma- 
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tions on the initial photon state. We suggest this as a 
tool to manipulate the photon distribution and make a 
filter in the state space of the radiation. From the litera- 
ture, however, we know that two-level systems with time 
dependent coefficients have played a central role in the 
understanding of quantum phenomena. For a review of 
such solvable time dependent models see the review [fj. 
Since a long time, many such models have been solved 
exactly and in this paper we are illustrating the filtering 
property by an application of these exact solutions. 
The Hamiltonian we consider is the following one 
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where the Pauli matrices are defined by 
[o-s,^] = ±2a ± 

0-3 I ±) = ± I ±). 



(2) 



The photon energy is given by Ml and the detuning is 
defined as 



Aw = LO — fi, 



(3) 



where htu is the energy separation in the two-level atom. 
The first term in the Hamiltonian (JTJ is a constant of 
the motion, and we can thus discuss the dynamics of the 
other terms neglecting this first one. As usual, its effect 
can be added later. Here we assume that the parameters 
are time dependent, git) and Aw(t). The latter is pre- 
sumably achieved by tuning the cavity frequency fl(t). 
We make an ansatz of the state in the form 

I*) = c o a_(0)|0,-) 

+ Er=i c « i a +( n ) \n-l,+) +a_(n) | n,->] . 

. (4) 

The initial states are assumed given by the coefficients 
c„ and aq_(n). After the interaction period, these latter 
ones go over into a±(n). 

Without the first term of the Hamiltonian QJ, we find 
for each family the equations of motion 



dt 



o+(n) 
a_(n) 
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a + (n) 
a_(n) 
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After a period of interaction, the solution of this problem 
can be written in the form 
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y^w^ e %Vn y/l — w n 1 [ a+(n) 

%Vn y/l — W n \/w^ J L a -( n ) 

(6) 

If we now, after the interaction, perform a measurement 
to find the two- level atom in one of the states | ±) , we 



find the cavity mode to be in the corresponding state 
with the photon distribution 



P+ = |a~(n + l)| 2 | C „ +1 | 2 
P- = |a-(n)| 2 | C „| 2 . 



(7) 



The asymmetry derives from our definition of the state 
, which allows the components in Eq. ijfjj to be labelled 
with the same n. With the atom emerging in the upper 
state, we need to start with one additional photon in the 
field. So the field distribution is modified by \a±(n)\ 2 , 
which we call filter functions. Before the observation of 
the state of the two-level atom, the normalization is given 
by 



1 = J2 P " = £ i P n + P n] = a_(0) H Co | 2 

n=0 n=0 
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[I a + (n) | 2 + |a-(n) | 2 



(8) 

After the measurement of the two-level state, the proba- 
bilities must be normalized to give 



(9) 



n=0 



For general initial conditions, the solution will de- 
pend on the phase tp n which can be quite complicated. 
However, the treatment simplifies if we restrict our at- 
tention to the two cases: 



Case (a) : \a°_ \ = 1; a% = 
Case(b): a°_ = 0, \a° + \ = 1. 



(10) 



In these situations, no interference terms need to be con- 
sidered. In the following the initial values are taken to 
be as in the case (a). Since the matrix in is uni- 
tairy it follows that the filter functions remain normal- 
ized, |a M (n)| 2 + |a^°(n)| 2 = 1. For the case (a), |a~(0)| 2 
does not change and hence the normalization in © be- 
comes J2n=0 l C »| 2 = L 

After the observation of the two-level state, we have 
projected out a modified photon state and applying a 
consecutive series of atoms with their initial states ac- 
cording to (|10(l . we can iterate the transformation (Q) 
without having to consider the phase variable. This al- 
lows us to enhance the filtering action exerted by a sin- 
gle atom affecting the cavity mode. It is important to 
point out the fact that each time an upper level atom 
is detected, and the photon distribution is modified by 
the corresponding filter function, the photon numbers are 
shifted by one unit as is seen in 0. 
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III. THE LANDAU-ZENER MODEL 

In the case solved by Landau and Zener, the parame- 
ters in Eq.JSJ) are given by 



.9 = .9o = constant 
Auj = 2Xt. 



(11) 



The solution integrated over the time (— oo, +oo) is given 
by 



w n = exp(— vn) 



Wo 
A ' 



(12) 



For the case (a) in (|10|) this gives 

| a + (n) | 2 = [1 — exp(— vn) 

| a„(n) | 2 = exp(— vn). 



(13) 



For the case (b) the roles of the two states are inter- 
changed. Thus for small values of the photon number, 
no transfer takes place, whereas for large photon num- 
bers, the atom approaches the adiabatic limit and the 
state is transferred. The filter action is based on this 
property. 

We can immediately see that the passage of a single 
two-level atom through the cavity effects a filtering which 
eliminates the low photon number states if the atom is 
found in its upper state and the high photon numbers 
if the lower state is observed. The experimentalist has, 
of course, no control over which is to be the case, but 
once the observation is done, the ensuing cavity state is 
determined. 

In order to observe the filtering action, we assume the 
cavity mode initially to be in a coherent state. Then we 
have 



exp( 



(14) 



From Eq.(7|), we find that projecting on the lower state, 
we have the normalized distribution 



exp(— ne ) : — 



(15) 

which is a Poisson distribution with the reduced average 
(n)_ = exp(— v)n . (16) 
Observing the other state, we find the distribution 



Pj 



c/v [1 — exp(— vn)] exp(— n)- 



(17) 



which is the difference between the original photon dis- 
tribution and the one with the reduced average (IS) , cn 
is the normalization constant. This is clearly a method 



to filter out the low photon numbers. The normalization 
coefficient can also easily be determined from the relation 



P + = J2 P n = ^{1 - exp [n(e-* - 1)] }. 



(18) 



After the detection of an upper level atom, the average 
photon number is 

/ \ 1 - e^p[n(e' v - 1) - v] _ 

\ n )+ = — r~, m — n—1. (19) 

w+ 1 - exp[n(e~ v - 1)] v ; 

If n 3> v the factor in front of fi is approximately unity 
and (n) w n — 1, while if v is large the factor becomes 
greater that unity and we may have (n) > fi. 

If we repeatedly observe the lower state after a se- 
quence of to atoms introduced in the lower level, we ob- 
tain the result 



(fie~ vm Y 

p-(m)=«p(-ne- wl )i ^ 



(20) 



which will greatly enhance the filtering action. The cor- 
responding distribution for the upper level is 



P+(to) cx 



f[\a™(n + v)\< 



cxp(— n) 



(n + to)! 



(21) 



As each observation is totally random, the achievement of 
a successful series of to projections will become smaller 
and smaller. However, the fact remains that after the 
observational sequence has been recorded, the state of 
the cavity mode is known. This holds independently of 
the actually observed sequence. If this contains alternat- 
ing upper and lower level observations, the appropriate 
transformation is to be taken from l|13|) . As the subse- 
quent projections have to be normalized, the normaliza- 
tion will contain the normalization coefficients of each 
preceding step in the process. As long as each step in- 
volves a projection on to one of the atomic states, no 
phase dependence needs to be taken into account. 

In order to illustrate the filtering behaviour of the 
Landau-Zener model, we show in figure 1, the filter func- 
tion | a°°(n) \ 2m for the repeated number of atoms, 
to = 1,5 and 25 and v — 0.126. For one atom we have 
the simple result, which does not provide an efficient fil- 
tering, but with increasing number of atoms, the high 
photon number parts become efficiently suppressed, thus 
pushing the average to lower and lower values. This may, 
for example, be used to cool an initially too hot thermal 
photon source. A projection on the upper level efficiently 
selects the high photon numbers instead. 



IV. THE DEMKOV-KUNIKE MODEL 

A. Results of the model 

In order to simplify the discussion, we consider only 
the case (a) of Ea. (|10|) . The parameters of the model 
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FIG. 1: This shows the filtering action achieved by the 
Landau-Zener model, when the outgoing atom is projected on 

to its lower level. The curves are presented for v = ^-f 2 -^ = 

0.126. The three curves give the result after the passage of 
one atom (m = 1), five atoms (m = 5) and twenty-five atoms 
(m = 25), which is indicated in the figure. We can easily see 
the increased efficiency of the filtering obtained. 



given in (jSJ are taken as 



Aw(t) = 2A t&xih(±) 
g(t) = gosech(^) . 



(22) 



This is a three parameter model, which integrated over 
the time (-co, +oo) gives the solution 



cos' 



cosli (ttTAo) 



o~(n) | 2 = - 
a°?{n) | 2 = 1- | a°?{n) \ 



(nTy/(g*n-A$) 



(23) 



For photon numbers such that g 2 n < Aq, we have 



cosh' 



(nTy/(At-g%n) 
cosh 2 (7rTA ) 



(24) 



as n — > 0. This must hold for all models as long as we 
consider the case (a) of Eq. I|1(J|I . The filter function pro- 
vided by this model, thus starts from unity, and within 
a range of (Ao/go) 2 of photon numbers goes to an os- 
cillating function with the amplitude sech 2 (7r7\4o). For 
small coupling constants go , this is thus essentially a low 
photon number filter like the Landau-Zener case. 

For a strong coupling case, go > Aq, the function is 
essentially an oscillating one, with the first maximum of 
\a™{n)\ 2 at 



n M 



+ 

T 2 g 2 



T 2 9 2 q 



(25) 



Even if the distance to the surrounding zeros is of the 
same order of magnitude, the trigonometric dependence 
on n provides a rather sharp cut-off around the maxi- 
mum. Application of the transformation (|24|l to a co- 
herent state with most probable photon number around 
riM will result in a considerable narrowing of the photon 
distribution. It is advantageous to choose the parameter 
Aq as small as possible, because then the oscillation am- 
plitude sech 2 (7rTAo) is maximized. The same happens 
when T is chosen short, however, then the oscillating pe- 
riod becomes large. We thus want to work close to the 
non-adiabatic limit, of the two-level model, in this case. 

In the non-adiabatic limit, the projection on the lower 
state of the atom gives the filtering 

| o~(n) | 2 = cos 2 (nTy/(g$n - A 2 ,)) sech 2 (7rTA ) 
« cos 2 (nTyfgJ^J , 



which gives a simple filtering around the maxima 
k 2 



n M 



T 2 gl 



(k = 0,1,2...) 



(26) 



(27) 



The zeros of the distribution are given by the expression 
«o = { ^M- (28) 



T 2 g 2 



o 



The width of the distribution is hence approximated by 

k 



An ~ n — tim 



(29) 



For a Poisson distribution, An ~ \fn~. Thus to set the 
maximum um to a large value requires a large k. The 
width (|29|l then scales as the square root of n with the 
implication that it matters little, which peak is chosen 
for the filtering. The factor Tgo can well be chosen close 
to unity in this parameter range. One may think that 
for optimal filtering action, the first maximum should be 
chosen. However, this does not exclude large values of 
Um - By a suitable choice of T, we can make both AqT and 
Tgo small even if go > Aq. Thus the three parameters 
may be utilised to achieve desired filtering operation. 

By sending atoms repeatedly through the cavity, in the 
case of a sequence of m recorded lower level observations, 
we obtain the filtering 



Pnim) 



,2m 




(30) 



The powers of the trigonometric function becomes 
rapidly sharper with increasing m. Thus the filtering is 
fast and efficient. As the repeated application of the 
filtering process sharpens up the photon distribution, 
we expect the states to rapidly become non-classical. 
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In order to explore this, we have investigated the Q- 
parameter, see 



Q 



(31) 



as a function of the coupling strength. For a range of 
initial coherent states, we find indeed negative values; 
see the next section. 

In the adiabatic limit, T growing large (i.e. AqT being 
large), the oscillational amplitude becomes small, and we 
expect the filtering process to become less efficient 



sech 2 (7rTA ) ~ 4exp(-27rTA ) < 1. 



(32) 



Then | a°°(n) | 2 from l|23|) decreases. In this case, we 
may also choose the weak coupling situation go < Aq 
and write 



cosh' 



cosh 2 (tt AqT) 



{exp [2nT (y/(A*-g*n) - A )] }™ | c n | s 
(-n^n) | c n | 2 . 



exp 



This form simplifies some of the arguments below. 



(33) 




FIG. 2: Fhis figure shows the filter functions for the Demkov- 
Kunike model. In our dimensionless units we have the param- 
eters, Aq — T = 0.1 and go = 4 . Here we also look at one, five 
and twenty-five atoms passing the cavity (m = 1, 5, 25). In (a) 
the lower level filter function |of°(n)| 2m is plotted. The in- 
crease of atom number rapidly narrows the distribution, and 
we can also see how the width scales approximately as the 
position of the peak of the function. This allows us to use the 
filter on a coherent state distribution situated at any peak of 
our choice. The filter function fT ! \ a, + {n + ^)| 2 is 

shown 

in (b). The maxima in (b) are shifted to the left by approx- 
imately m/2 units. The numbers in the figures are again 
indicating how many atoms m have been detected. 



B. Numerical results 

In this section we illustrate our argument by some nu- 
merical examples. As all relevant parameters has the di- 
mension of time or frequency, fixing one time scale allows 
us to use dimensionless units in the calculations. 

As we have seen, if we repeatedly project on the lower 
or the upper atomic level, the new photon distribution 
will be modified by \a+{n)\ 2m or Jl^Li \ a +( n + v )\ 2 rc ~ 
spectively. In the latter case the initial photon distribu- 
tion will also be shifted by m units. In figure 2 (a) and 
(b) we plot the two functions above for m — 1,5 and 25. 
The maxima in (b) are shifted to the left by m/2 units. 
For small AqT the function |af°(n)| 2m will be approxi- 
mately zero everywhere except where n « tim and tim is 
close to the maxima (12711 . If the initial photon distribu- 
tion, given by \c n \ , has a maximum near one riu and a 
width smaller than the period of the filter function, then 
we expect a sharp final distribution if a series of m lower 
level atoms is detected. We know that the Q-parameter 
defined in l|31[l is positive for a super-Poissonian and neg- 
ative for sub-Possonian state, which signals non-classical 
light. A plot of Q as a function of go is shown in figure 
3 for an initial coherent state with n = 100, the number 
of atoms m = 25 and T = Aq = 0.1. We see that the 
field distribution P~ is sub-Poissonian for a large range 
of gQis. Note that one could use the Q-parameter to find 
the optimal parameters T and go in order to obtain the 



most non-classical light. As noted in the previous sec- 
tion, since the width of the maxima of the filter function 
goes as yJriMi just like the width of a coherent state, the 
ratio between them is independent of n. So if we know 
go , T and Aq , the filtering effect is chiefly independent of 
which maximum um the initial coherent state is centred 
at. The distribution P~ for an initially coherent state, 
centred around the third maximum n — 25 of the filter 
function, is shown in figure 4 with T = Aq = 0.1, go = 4 
and the number of atoms m = 1, 5 and 25. 

In order to have a low photon number filter it was 
found to be advantageous to use a large adiabaticity pa- 
rameter A = AqT and a small coupling constant go < Aq. 
A large adiabaticity parameter means that the oscillating 
part of the filter function is small and the weak coupling 
implies that the hyperbolic behaviour in the filter func- 
tion extends over a long distance. In figure 5 the filter 
function |a^°(n)| 2 is plotted for the parameters go = 0.2, 
T = 0.9 and Aq = 0.6. Except for the oscillations, this 
should be compared with the case m = 1 in figure 1. 



C. Analytical approximations 

We have considered two limits of the two-level model: 
The non-adiabatic one which could give a sharpened pho- 
ton distribution and the adiabatic one which is used to 
filter out the high or the low photon number states. In 
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FIG. 3: This shows the variation of the Q-parameter with 
coupling strength. We must have Q > — 1, but all negative 
values indicate non-classical photon statistics. The intial state 
assumed here has average photon number n = 100; the inter- 
action is defined by T = Ao = 0.1. In this figure, we consider 
the case with m = 25 atoms having successfully passed the 
cavity and been found in their lower level. The insert zooms 
in on the first non-classical regions, and shows the degree of 
non-classicality which is achieved. The limit for classical be- 
haviour, Q — , is indicated by the dotted line. 




FIG. 4: This figure shows the narrowing of the phton distri- 
bution actually achieved with one, five or twenty-five atoms 
(m = 1, 5, 25) acting on a coherent state centered on n = 25, 
which in this case is the second maximum; c.f. figure 2 (a). 
The dimensionless parameters are T = Ao = 0.1 and go — 4. 



this section we derive analytic expressions for the widths 
of the sharpened distribution and the low photon number 
filter. 

The filter function is given by Eq. (|26|l in the non- 
adiabatic limit, and its maximum by the result (|27|l . Af- 
ter the passage of m atoms, all of them detected in their 
lower level, we consider the photon distribution. The 



FIG. 5: This is the filtering function after the passage of 
one atom in the limit when the Demkov-Kunike model acts 
as a low photon number filter. Here T = 0.9, Ao = 0.6 and 
go = 0.2. Except for the oscillations, this is similar to the 
Landau-Zener result, figure 1 (m = 1). 



FWHM, An Q , of the filtering function is determined from 
cos 2 " 1 (nTgoy/n M ± An a /2) = 
cos 2m (n M ± arccos (O^ 1 / 2 ™))) , 



(34) 



where the signs ± derive from the fact that the filtering 
function is not symmetric around tim due to the y/n - 
behavior. Equation (|34|) determines two different An a 
and to get the appropiate value we take their average 



An n 



(35) 



where x m = arccos ( 2 \/0.5) which starts at 7r/4 and falls 
to zero for increasing m. For an initial Poission distribu- 
tion, we have for the half-width to a good approximation 
An p — \JSum m (2), and the width of the final distribu- 
tion P~ (m) is approximately 



An 2 = 



Wx 2 m n M ln(2) 

An 2 + An 2 " 7r 2 T 2 .g 2 ln(2) + 2x 2 



(36) 



For a narrow An a we have Art An a ~ x m . A plot 
of An for g = 4.5, A = T = 0.1 and n = 19.8 as a 
function of m is shown in figure 6 together with exact 
numerical results for the FWHM. 

Now we turn to the adiabatic limit and the low photon 
number filter. With the right parameters it is possible 
to achieve such filtering by repeatedly projecting on the 
lower atomic level. The width, An at e _1 , of this low 
photon number filter is derived from 



cosh 2 " 1 irTJA 2 - 5 2 An sech 2m (7rTA ) = 



(37) 
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FIG. 6: This figure compares the approximate expression 
l|H6> (solid line) for the width of the final photon distribution 
with the exact one obtained numerically (diamonds). The 
dimensionless parameters are here T = Aq = 0.1 and go — 4.5. 
The plot shows the influence of increasing atom number m on 
an initially coherent state with n = 19.8. 



which determines An. The result is complicated and not 
very informative. However if we instead use the fact that 
we are in the adiabatic regime and take the last expres- 
sion in Eq. (|33[1 as the filter function, we obtain the 
width 



An 



A 



(38) 



The width clearly decreases as m becomes large, so that 
the high photon numbers are effectively filtered out. In 
Fig. 7 the width (|38ll is plotted together with the exact 
width obtained from Eq. I|37|) as a function of m for the 
parameters Ao = 2 and T — go = 1. Again we have a 
very good agreement. 



V. CONCLUSION 

In this paper, we consider the situation where individ- 
ual atoms are made to traverse a cavity at the will of the 
experimenter. During the passage of the atom through 
the cavity, the parameters of the model are made to vary 
in order to achieve a desired effect. This variation must 
be slow enough to preserve the energy levels of the inter- 
acting systems intact. For slow atoms this should present 
no problem, but the time intervals involved are limited to 
those allowed by the various relaxation processes which 
are omitted in our treatment. When actual experimen- 
tal parameters are introduced, the influence of relaxation 
processes can easily be included by performing numeri- 
cal integration of the time evolution. Here we want to 
present the main ideas, and such detailed investigations 
have been excluded. 




FIG. 7: This figure compares the approximate expression I38H 
(solid) for the width of the final photon distribution with the 
exact one obtained from 13711 (diamonds). Here we are in 
the limit where the interaction acts as a low photon filter 
and we see the influence of increasing atom number m. The 
dimensionless parameters are Aa = 2 and T = go = 1. As we 
have a low photon filter, this may be interpreted as the range 
over which the filter acts. Alternatively we may use it as a 
high photon one with the corresponding interpretation. 



In order to illustrate our ideas, we choose to utilise the 
existence of analytically solvable models for the interac- 
tions. Choosing the simple Landau-Zener model as a pro- 
totype, we then investigate in detail the behaviour for the 
related but much more realistic Demkov-Kunike model. 
They both allow filtering action for either low or high 
photon numbers. If low photon numbers are favoured, 
the process can be used to cool the photon distribution, 
but favouring high photon numbers we can sharpen the 
distribution or heat it up. For sharpening it, the Demkov- 
Kunike model provides a much more efficient tool, which 
we investigate in some detail. In particular, we show that 
strongly sub-Poissonian distributions can be achieved. 

There are many treatments of time-dependent two- 
level systems in the literature. Many cases have been 
solved analytically and the standard method is to reduce 
the problem to a differential equation with a known solu- 
tion ■ Such models have been utilized widely to model 
time-dependent atomic processes Any one of these 
could be used in our present scheme to analyze the re- 
sult of the atom-cavity interaction. However, as we are 
here interested in the modification of the photon distri- 
bution, these earlier works are only of interest so far as 
they provide novel filter functions for the quantum field. 
Consequently we have here choosen to investigate two 
of the best known models only; they provide enough in- 
sight into the possibilties and limitations of our suggested 
approach. More intricate technical details of the actual 
solution are of little interest in the present context. 

In order to obtain the most dramatic effects, we have to 
let several atoms pass through the cavity. Their state is 
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determined experimentally after the interaction period, 
and the corresponding modification of the photon dis- 
tribution becomes known. This occurs for either of the 
two possible outcomes of the experiment. However, if we 
want to achieve a pre-assigned series of observations, this 
can only be found with a certain probability. Several at- 
tempts may be necessary to observe the desired sequence. 
However, once it has been found, the cavity state is also 
known, and it may be used in any subsequent experiment 
within the period allowed by the cavity relaxation rates. 

Our suggested approach is, of course, part of the broad 
range of works called " state preparation" and " state en- 
gineering" . Much work has been done in this area, and 
we cannot possibly relate our effort to all of these. For 
an overview consult e.g. the special issue |2j on Quantum 
State Preparation and Measurement. Here we only want 
to point out, that the modification of the state is de- 
termined by the dynamic interaction between the atoms 
and the cavity mode. The repeated observation leading 
to projection onto the new cavity state does not, in it- 
self, shape the distribution, it only picks one out of two 
possibilities. 

There are many cavity QED schemes where the field 
in the cavity is modified by transit of atoms through the 



cavity, see e.g. P3 and H2 and refereces therein. Meth- 
ods for prepering various states of the cavity field have 
been proposed, e.g. Fock states Schrodinger cat 

states [K| and EPR states between two cavity modes 
|l3j. In these works, the time dependent feature is the 
transit of the atom through the cavity, and a suitable 
measurement at the end of the interaction probes the re- 
sult achieved. The idea is, of course, closely related to 
our approach, but we have added the novel feature that 
the parameters of the coupled atom-cavity system are 
deliberately changed so that we achieve a desired effect. 
Thus the modification of the quantum state is effected by 
the control imposed by the experimenter and not only by 
the interaction followed by a detection procedeure. 

We have utilized solvable models to illustrate out ba- 
sic ideas. These make it possible to assert the influence 
of the various parameters in a straightforward way. In 
an actual experimental situation, we may, of course, per- 
form the corresponding dynamical calculations on a com- 
puter, which allows one to introduce various effects re- 
lated to experimental conditions and non-ideal observa- 
tional methods. Such considerations are, however, best 
performed in connection with realistic laboratory situa- 
tions. 
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